Above threshold s-wave resonances illustrated by the l/2 + states in 9 Be and 9 B 
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We solve the persistent problem of the structure of the lowest l/2 + resonance in 9 Be which 
is important to bridge the A = 8 gap in nucleosynthesis in stars. We show that the state is a 
genuine three-body resonance even though it decays entirely into neutron- 8 Be relative s-waves. The 
necessary barrier is created by "dynamical" evolution of the wave function as the short-distance 
a- 5 He structure is changed into the large-distance n- 8 Be structure. This decay mechanism leads 
to a width about two times smaller than table values. The previous interpretations as a virtual 
state or a two-body resonance are incorrect. The isobaric analog l/2 + state in 9 B is found to have 
ener 
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energy and width in the vicinity of 2.0 MeV and 1.5 MeV, respectively. We also predict another 
9 B with similar energy and width. 



PACS numbers: 21.45.-v, 21.60.Gx, 27.20.+n 
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Introduction Bound states, resonances and other con- 
tinuum states are well understood and described for two 
particles interacting through a potential. By weakening 
the attraction of the potential bound states are pushed 
upwards into the continuum as resonances or virtual 
states when a barrier is present or absent, respectively 
[![. For neutral particles virtual states arise for s-waves 
whereas resonances emerge for higher partial waves. De- 
creasing the attraction further until the resonance energy 
is above the potential barrier leads to an increase of the 
resonance width. Correspondingly the related ^-matrix 
pole moves in the complex energy plane as a resonance 
with non- vanishing imaginary part. 

For three particles interacting via two- and three-body 
potentials the continuum structures can be much more 
complicated due to combinations of the different struc- 
tures for the three two-body subsystems 0, 0| • One in- 
triguing possibility arises when one of the two-body sub- 
systems has a low-lying s-wave resonance produced by a 
confining Coulomb barrier, and the third neutral parti- 
cle has dominating s-wave attractions from the first two 
particles. Even when all higher partial waves are vanish- 
ingly small, the structure of the three-body continuum 
state is a priori not easily determined or described. 

Let us assume that the two-body resonance is very nar- 
row (long-lived) and the three-body energy is above zero 
but less than the two-body resonance energy. Then the 
three-body continuum state resembles a two-body bound 
state of the third particle and the composite resonance of 
the two first particles. The lifetime would be determined 
by the lifetime of the two-body resonance. When the 
three-body energy is pushed upwards above the two-body 
threshold by decreasing the attraction, the corresponding 
structure can be described as a two or three-body virtual 
state or resonance [H, H| . 

The purpose of the present letter is first to determine 



in general which structure arises, and second specifically 
to solve the long-standing controversy of the 9 Be(l/2 + ) 
continuum state. This state is important in the nuclear 
synthesis of light nuclei in stars [4-6], and has therefore 
received lots of attention both theoretically (Bi^ and ex- 
perimentally [mHUI • In a measurement of photo disinte- 
gration the cross section is interpreted and parametrized 
via i?-matrix analysis as one neutron and the 8 Be ground 
state in a two-body s-wave resonance [llfl . This seems to 
be against the two-body quantum mechanical description 
as such a state cannot survive as a resonance. In another 
interpretation the same neutron- 8 Be system is described 
as a virtual state [§] but the resulting cross section does 
not reproduce the measurement [ill ]. 

The 9 Be(l/2 + ) structure is most often assumed to be 
one neutron and the 8 Be ground state [l(| but some- 
times also the a + a + n recombination reaction is as- 
sumed to proceed by a-capture on the 5 He ground state 
[0] . It is apparently very difficult to avoid assumptions of 
two-body sequential structures and processes via subsys- 
tems of either 8 Be or 5 He. Interestingly a two-center 
Born-Oppenheimer model based on symmetries alone 
may combine these structures as in where the 1/2+ 
is lowest at large distance whereas a 3/2~ state is low- 
est at small distance. We shall allow an entirely general 
three-body structure without a priori assumptions of sub- 
structures or decay mechanisms. 

Formulation Let us consider three composite struc- 
tures as point-like particles denoted n, ai and ct2- The 
two mass scaled Jacobi vector coordinates, (x, y), can be 
substituted by hyperspherical coordinates {p, a, ft x , 
where (Q x ,fl y ) describe the directions of (x,y), p = 
\/x 2 + y 2 and a — arctan(x/y), see 0]. We solve this 
three-body problem by use of adiabatic hyperspherical 
expansion of the Faddeev equations, i.e. the angular- 
equations are solved for each p, providing a set of an- 



2 



gular eigenfunctions <p n and their corresponding eigen- 
values A n . The three body wave function is then written 
as ip = -F72 J2 n fn(p)4'n(x, y), where n labels each of the 
adiabatic terms. The radial functions f n (p) are obtained 
after solving a coupled set of radial equations where the 
X n angular eigenvalues enter as effective adiabatic poten- 
tials. The coupling between the different adiabatic terms 
appears through the functions P nn ' (p) and Q nn ' (p) de- 
fined for instance in |3|]. 

Each adiabatic potential describes a specific relative 
structure of the three particles for a given root mean 
square radius, p, with wave function <f> and eigenvalue 
A. When only one adiabatic potential is considered, the 
coupled set of radial equations reduces to 

cP \{p) + 15/4 2m(E-V 3b (p)) ] 

- W 2 + — ? Qip) p J f(p) = 

where Q is the diagonal coupling term Q nn , which is in 
general not zero (contrary to what happens with P nn i, 
whose diagonal terms are zero). V^bip) is the three-body 
potential usually used in three-body calculations to take 
into account all those effects that go beyond the two- 
body interactions. The total wave function, ip, and Q 
are given by 

V> = -^f(p)4>(x,y) , Q(p) = {<t>\^\4>)n ■ (2) 

The expectation value is over angular coordinates, f2, ex- 
cluding only p. When only s-waves contribute for both 
x and y the angular wave function only depends on p 
and a. For short-range attractive two-body interactions 
the angular eigenvalue \{p) would be monotonously in- 
creasing towards a constant asymptotic value. 

Narrow two-body resonance. When the ot\ — a 2 two- 
body s-wave interaction supports a bound state, the adi- 
abatic potential approaches the bound state energy at 
large distance. This reflects a two-body structure corre- 
sponding to particle n far away from the a\ — a 2 bound 
state. For less a\ — a 2 attraction this state moves into 
the continuum. With a confining Coulomb barrier from 
repelling charges on the particles, the state would appear 
as a resonance at low energy E 2 . The adiabatic poten- 
tial would then approach the positive value equal to E 2 
corresponding to the large-distance two-body structure 
of n far away from the resonance a% — a 2 . We illustrate 
in FigO] by the specific examples, 9 Be(a + a + n) and 
9 B(a + a +p). 

This description is only correct when the two-body res- 
onance width T 2 is very small and the coupling to the 
three-body continuum states can be ignored. In general 
the first adiabatic potential is crossed by numerous oth- 
ers while p increases. However, the couplings are neg- 
ligibly small for three-body energies E until distances 
p ~ 9y/E - E 2 /T 2 where the energies are in MeV and p 
in fm. Thus for small T 2l say eV or keV, the couplings 



for moderate energies below f MeV can be neglected far 
outside the distance where the short-range n — [a\ — a 2 ) 
interaction has vanished. Then the system can effectively 
be described as a two-body system until the two-body 
resonance eventually decays. 

Let us now consider the three-body system with inter- 
actions leading to a three-body state of positive energy 
E r but below E 2 . Effectively this is a bound (ai — a 2 ) 
two-body state, or rather a resonance decaying precisely 
with the width T r — T 2 of the u\ — a 2 resonance. For 
less n — (ai — a 2 ) s-wave attraction this two-body bound 
state moves into the continuum above E 2 . The expecta- 
tion is that the proper description is as a virtual state 
with no width in contrast to a resonance However, 
this is not necessarily true. 
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FIG. 1: The two lowest energy levels for 9 Be (a), and 9 B 
(b), and the resonance energies of the corresponding two-body 
subsystems [T3|. For 9 B the quoted f/2 + state corresponds 
to the estimation obtained in this work. The widths of the 
resonances are represented by the shadowed regions. 

The 9 Be(l/2 + ) example. The low-lying states of 9 Be 
are expected to be well described as cluster states consist- 
ing of one neutron and two a-particles [l6| ■ The a — a in- 
teraction, including short range attraction and Coulomb 
repulsion, produces a low-lying s-wave two-body reso- 
nance at 0.0918 MeV with a width around 9 eV. Adding 
one neutron in s-waves leads to angular momentum and 
parity 1/2+ of the resulting three-body system. Such a 
state is listed in the tables of energies [17[ at an excita- 
tion energy of 1.684 MeV, or 0.110 MeV above threshold, 
with a width of 0.217 MeV (see the upper part of FigQJ,). 
Thus the state is above the two-body resonance energy 
by 0.018 MeV and s-waves are most likely the dominating 
composition. 

The listed width is much larger than the distance to the 
two-body resonance threshold and even about two times 
larger than the three-body energy itself. It is a peculiar 
resonance structure which apparently extends into the 
bound state region below the threshold. These values are 
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consistent with photodissociation cross section measure- 
ments and the entangled i?-matrix analysis of an a priory 
assumed resonance [U[ . The fitting parameters are posi- 
tion and energy dependent width of a two-body neutron- 
Be resonance of s-wave character. Several years prior 
to this analysis it was suggested that the state should be 
understood as a virtual state and the photodissociation 
cross section correspondingly analyzed [9(. However, this 
does not reproduce the measurements in [ill ]. 
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An attractive pocket and a constant large distance po- 
tential without any barrier in the transition region is not 
able to support a resonance of finite width at energies 
above the large distance asymptotic value. However, in- 
clusion of the coupling term Q (see Eq.(JTJ)), provides an 
otherwise totally absent barrier as seen in the inset of 
FiglU This all decisive barrier then arises from a strong 
p dependence of the intrinsic (angular) wave function </>, 
as seen from the definition in Eq.@. 
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FIG. 2: Lowest adiabatic potentials for 9 Be and 9 B as a func- 
tion of the hyperradius. The inset shows the 9 Be lowest po- 
tential with and without the rearrangement coupling term Q. 

Three-body results for 9 Be(l/2 + ). The shortcomings 
of many previous methods are the initial assumption of 
two-body character of this structure. Since the final de- 
cay products are three particles we proceed to treat the 
system as a three-body system. We use the well estab- 
lished nucleon-nucleon and a-nucleon interactions from 
11,0. 

The adiabatic potentials produced through the adia- 
batic hyperspherical expansion method for these quan- 
tum numbers are shown in Fig|2] The higher ones show 
pronounced peaks at about 15-18 fm. The origin of the 
peaks is in the crossing between different angular eigen- 
values. In particular, the Q-functions in the effective 
potentials, (see Eqs.dJ) and flU)), are responsible for the 
behaviour shown in the figure in the vicinity of the cross- 
ings. These couplings involve second derivatives of the 
adiabatic eigenf unctions, and they therefore reflect re- 
structuring of these functions. The lowest potential has 
a dominating attractive pocket at small distance. After 
hyperradii p larger than 12 — 14 fm the potential stabi- 
lizes at the resonance energy of 0.091 keV for the 8 Be 
ground state. This stable region continues until inter- 
rupted by the crossings with the (infinitely many) higher 
potentials. The first of these crossings occurs at about 
p = 130 fm. 



FIG. 3: The partial wave decomposition of the lowest adia- 
batic angular wave function for 9 Be (thick) and 9 B (thin) as 
function of hyperradius p. The partial angular momenta l x 
and l y correspond to the coordinates indicated in the figure. 



For l x = l y = 2 and l x = 
can not be distinguished. 
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3 the curves for Be and B 



The three-body restructuring can be seen from FigJ3] 
where the small distance structure of p less than 9 fm is 
p- waves between neutron and ev-particles. This is due to 
the p3/2-resonance which provides the main part of the 
attraction. As p increases above 10 fm this partial wave 
is rapidly substituted by s-waves which are energetically 
more advantageous at larger distance where the attrac- 
tion vanishes and the centrifugal barrier dominates. At 
much larger distances an increasing number of partial 
waves contribute corresponding to the neutron far away 
from the two-body system of two a-particles in the spa- 
tially much smaller s-wave resonance, i.e. 8 Be in the 
ground state. 

The combined result of the lowest adiabatic poten- 
tial and the diagonal coupling is able to support a res- 
onance. We compute the energy and width numeri- 
cally as the S'-matrix pole found by complex scaling 
20]. By adding a structureless short-range potential 
(Vsb = Vq cxp(— p 2 /pq), with pq = 5 fm), which con- 
tributes only in the pocket region, we can move the res- 
onance energy by modifying the Vq strength but without 
disturbing the resonance structure. This is useful both 
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because the three-body computation can not place the 
resonance at precisely the correct measured energy, and 
because the width is strongly dependent on the height 
and thickness of the barrier at the correct energy. 
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FIG. 4: Width of the resonances for 9 Be and 9 B as function 
of the energy which is varied through the strength Vb of the 
three-body potential. The solid and dashed curves are the 
WKB results with a knocking rate corresponding to To = 
0.6 MeV (r = T e~ 2S ) for both nuclei The dot-dashed 
curve results from complex scaling for 9 B. The square and the 
circle are from the R-matrix analysis in [ll[ and the table in 
[T3|, respectively. The down triangle is obtained by direct fit 
of the cross section in The triangles at about 2 MeV are 
the first and second resonances of 9 B. 

In FigH] we show the width as a function of the res- 
onance energy. The solid line shows the WKB estimate 
for 9 Be. For an energy just above the threshold energy of 
the 8 Be ground state (0.0918 MeV) the width is about 0.1 
MeV. This agrees with the result obtained by fitting the 
measured peak in the photodissociation cross section 11 1 
(triangle down). The width increases slowly up to about 
0.7 MeV at the top of the barrier. At smaller energies 
the width is vanishingly small due to the thick barrier 
provided by the 8 Be structure. 

The complex scaling computations for 9 Be present nu- 
merical difficulties due to the fast increase of the width 
for energies above 0.0918 MeV. To get an estimate we 
have instead attached a unit charge to the neutron, which 
immediately leads to the results shown in Fig(4] for the 
l/2 + analog state in 9 B. A continuous decrease of the 
charge leads us back to 9 Be which is approached and 
finally reached by extrapolation. Unfortunately the ac- 
curacy only allows the conclusion that the width at an 
energy of 0.11 MeV is larger than 0.06 MeV and most 
likely around 0.1 MeV but 0.2 MeV is not numerically 
excluded. The inaccuracy is due to the very large width 
compared to the distance to the threshold. Then the 
poles computed by the complex scaling method are very 



difficult to distinguish from the continuum background. 
This becomes increasingly worse as the charge of the neu- 
tron is decreased from unity to zero where our present 
techniques prohibits a clean result. 

Comparison to widths obtained in previous works re- 
quires precision in the definitions of a resonance. The 
ambiguities arise when the resonance is broad and then 
necessarily asymmetric. To account for the energy de- 
pendence of the decay probability the i?-matrix theory 
employs an energy dependent width. This immediately 
implies that any number claimed to describe the width 
must be an average of some kind. The "observed" width 
in i?-matrix theory is in [22] most directly related to the 
full width at half maximum, or the S'-matrix pole, and 
then adopted as the width in tables of resonance proper- 
ties [I3|- The results quoted in 11 1 and are shown 
in FigfJlby the square and the circle, respectively. 

Our computed three-body resonance width is from the 
WKB tunneling and crude extrapolations from the imag- 
inary value of the pole of the S-matrix. Our estimate is 
very inaccurate but we expect a value of about 0.1 MeV 
which is only about half the "observed" table value from 
.R-matrix theory. This rather large discrepancy can be 
due to inaccuracies in the complex scaling extrapolation 
to zero charge of the neutron, or the WKB approxima- 
tion of only one adiabatic potential combined with the 
uncertainty in the knocking rate estimate. However, we 
believe that the main reason is that our three-body bar- 
rier, which entirely is responsible for the width, arises 
from three-body restructuring effects inherently impossi- 
ble to include in the i?-matrix analysis. An unambiguous 
settlement of this accuracy issue would require the use of 
another method dedicated to precise width computations 
near threshold. 

Attempts to settle the issue experimentally face the 
problems inherent to relatively very broad resonances. 
Population in a reaction or beta-decay provide informa- 
tion about lifetime, which necessarily is an average, or de- 
tails about decay probability as function of energy, which 
requires a model for analyzing the data. The direct mea- 
surements in photo dissociation [ill ] reveal an asymmet- 
ric peak with a width of about 0.1 MeV. This is con- 
sistent with corresponding computations in the present 
model with a method to compute strength functions by 
discretization of the three-body continuum • 

Three-body results for 9 B(l/2 + ). The isobaric analog 
state should exist in 9 B but so far it has never been found. 
The influence of the additional Coulomb interaction in 
the present fragile case could be substantial and perhaps 
even destructive. The corresponding low energy three- 
body and two-body levels are shown in FigfTJ). To discuss 
this problem we also show in Fig|2]the corresponding low- 
est adiabatic potential for 9 B. The additional Coulomb 
repulsion due to the proton in 9 B has a relatively small 
effect. At small distances the pocket is as pronounced but 
almost entirely above zero energy. The barrier is higher 
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and thicker at small energy where the adiabatic potential 
itself already exhibits a small barrier. The same constant 
energy is approached at larger energies again correspond- 
ing to 8 Be in addition to the Coulomb tail from the pro- 
ton. The partial wave decomposition in FigJ3] resembles 
the 9 Be results with a tendency towards an increase of 
higher partial waves due to the additional Coulomb po- 
tential. 

In Fig|4]we show the WKB estimate (dashed line) and 
the resonance width obtained after a complex scaling cal- 
culation (dot-dashed line). In both cases the width is 
small even above the 8 Be ground state energy due to the 
additional adiabatic barrier and the extended Coulomb 
tail. However, the precise value of the energy depends 
on the three-body potential used in Eq.JT]). A minimum 
value for the attractive strength of the three-body force 
can be obtained by placing the 9 Be resonance at the 8 Be 
resonance energy threshold. Using this three-body force 
we can then estimate a lower limit for the 1 /2 + resonance 
in 9 B, which is found to be slightly below 2 MeV with a 
width of 1.3 MeV (left triangle up in the figure). 

This resonance is accompanied by a second state at 
2.05 MeV with a width of 1.6 MeV (right triangle up in 
the figure). This second resonance is very stable, basi- 
cally independent of the structureless three-body force. 
Therefore, further decrease of the three-body attraction 
in V 3 b would eventually make this second resonance at 
2.05 MeV the first 1/2+ excited state. The two reso- 
nances are dominated by s and p-waves between proton 
and a-particle, respectively. Crudely speaking these two 
structures correspond to 8 Be+proton and 5 Li+a. Since 
the large-distance p-wave properties essentially are unaf- 
fected by the three-body potential the second resonance 
remains close to the same energy. Thus the unobserved 
resonance in 9 B may in fact be either two or a combina- 
tion of two resonances. In any case we expect genuine 
three-body structures with an energy around 2 MeV and 
a width in the vicinity of 1.5 MeV. 

Discussion and Conclusions. The l/2 + continuum 
structures of 9 Be and 9 B are computed as genuine three- 
body resonances. The energies are above the threshold 
for forming the ground state 8 Be-resonance and the par- 
tial wave decompositions are dominated by s- waves in the 
Jacobi coordinates connecting the two a-particles and 
their center-of-mass and the nucleon. Nevertheless the 
Faddeev component corresponding to the other Jacobi 
coordinate presents a very different partial wave decom- 
position where the a-nucleon p 3 / 2 attraction maintain p- 
waves at hyperradii smaller than 9 fm rapidly changing 
into s-waves at 10 fm. This "dynamic evolution" from 
a— 5 He at small distance to neutron- 8 Be at intermediate 
and large distance reconciles the two limits for the reso- 
nance structure. The restructuring of the wave function 
results in a potential barrier similar to an above barrier 
reflection at a discontinuity. 

The structures are then genuine three-body resonances 



mixing s and p-waves at small distances while at large 
distances turning into a two-body system entirely of s- 
waves for a nucleon and 8 Be in the unbound ground state. 
The deceiving appearance as a two-body virtual state 
is incorrect. The interpretation as a two-body nucleon- 
8 Be resonance is also incorrect since the small-distance 
structure is of genuine three-body character and this is 
the very reason for the existence of a barrier allowing 
the appearance of resonance features. Furthermore the 
parameters from i?-matrix analysis of the experimental 
data is misleading because of the incorrect but crucial 
assumption of the existence of a two-body nucleon- 8 Be 
resonance. 

The astrophysical naa recombination rate is proba- 
bly unaffected provided the corresponding cross section 
is obtained by precisely the same parametrization for the 
measured inverse process of photodisintegration. The 
problem only seems to arise when a different procedure 
is applied in these mutually inverse processes. However, 
all resonance decays do not necessarily proceed through 
two-body channels. 

The implication is that the decay mechanism is entirely 
through the 8 Be ground state as assumed in most previ- 
ous publications. However, the present tabulated value of 
the resonance width emerges through an averaging pro- 
cedure of data analysis and parametrization of a decaying 
two-body structure. The energy dependent width arising 
from the i?-matrix interpretation as a two-body structure 
is suspicious since the three-body effects responsible for 
the barrier and the width are not included in the analy- 
sis. The true resonance lifetime should be related to the 
tunneling probability through the barrier arising from re- 
structuring the three-body wave function. The width is 
more likely directly found from the peak in the measured 
cross section. We conclude that the controversy over the 
1/2+ continuum structures of 9 Be and 9 B is resolved in 
a full three-body model. 
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